Quantum phase transition in an atomic Bose gas with a Feshbach resonance 
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We show that in an atomic Bose gas near a Feshbach resonance a quantum phase transition occurs 
between a phase with only a molecular Bose-Einstein condensate and a phase with both an atomic 
and a molecular Bose-Einstein condensate. We show that the transition is characterized by an Ising 
order parameter. We also determine the phase diagram of the gas as a function of magnetic field and 
temperature: the quantum critical point extends into a line of finite temperature Ising transitions. 



Introduction. — One of the most important recent de- 
velopments in the field of uhracold atomic gases is the 
apphcation of Feshbach resonances. A Feshbach reso- 
nance in the scattering amplitude of two atoms occurs 
when the total energy of the atoms is close to the en- 
ergy of a molecular state that is weakly coupled to the 
atomic continuum In the alkali gases of interest, 

this coupling is provided by the exchange interaction, and 
consequently the magnetic moments of the two atoms 
and the molecule differ substantially. As a result the en- 
ergy difference between the molecule and the threshold of 
the two-atom continuum, known as the detuning 5, can 
be experimentally tuned by means of a magnetic field. 
Moreover, by sweeping the magnetic field from positive 
to negative detuning through the Feshbach resonance, it 
is actually possible to form molecules in the atomic gas 

iiiiQ. 

Very recently, it has even been possible to create 
a Bose-Einstein condensate (BEG) of molecules in an 
atomic Fermi gas with a Feshbach resonance 0, 0, ■ 
This achievement is of particular importance because it 
offers the opportunity for observing a Bardeen-Cooper- 
Schrieffer (BCS) transition in a dilute atomic gas [llLfl^ . 
The idea here is to start from a pure molecular conden- 
sate and then sweeping the magnetic field back to positive 
detuning to create a much colder atomic Fermi gas than 
could have been achieved without the intermediate step 
of forming a molecular Bose-Einstein condensate. It is 
important to realize that this idea strongly hinges on the 
fact that a smooth BEC-BCS crossover exists in this sys- 
tem as one changes 5 0, J^] . In this Letter we show 
that analogous experiments varying 5 in an atomic Bose 
lead to a true phase transition rather than a crossover 
(lif. The reason for, and the nature of, this phase tran- 
sition can be understood as follows. 

The argument contains two important ingredi- 
ents. The first ingredient is that the coupling be- 
tween the atoms and molecules in the gas is pro- 
vided by an interaction energy that is proportional 

to /dx (V'Ji(x)V'a(x)?Aa(x) -h'!/'l(x)^I(x)'!/'m(x)), whcrC 

'(/'a(x) and '(/'in(x) annihilate an atom and a molecule at 
position X, respectively [13, 0, 0|. Such a couphng 
implies that if the gas contains an atomic Bose-Einstein 



condensate, and therefore has a nonzero value of ('0a (x)), 
the gas must necessarily also contain a molecular Bose- 
Einstein condensate and have a nonzero value of (?/'m(x)). 
However, the reverse is not true, and it is possible for the 
gas to contain only a molecular Bose-Einstein conden- 
sate. 

The second ingredient is that the symmetries of these 
two different phases of the gas are different. In the nor- 
mal phase the gas is invariant under the phase trans- 
formations ■0a (x) e*^?Aa(x) and ?Am(x) e^'^'0m(x). 
The additional factor of two in the transformation of 
the molecular annihilation operator follows also from 
the above interaction energy and is physically related 
to the fact that an molecule consist of two atoms. If 
the gas contains both an atomic and a molecular Bose- 
Einstein condensate (AC-t-MC) this U{1) symmetry is 
completely broken, and no residual symmetry exists. 
However, if the gas contains only a molecular conden- 
sate (MC), a residual discrete symmetry remains because 
('0m(x)) — > ('i/;m(x)) for 6 = TT. This phase therefore only 
breaks the U{l)/Z2 symmetry spontaneously. As a re- 
sult there must exist an Ising-like phase transition be- 
tween the state with only a MC and the state with an 
AC-I-MC, in which the residual Z2 symmetry is sponta- 
neously broken. Note that in the case of an atomic Fermi 
gas the BCS phase is characterized by a nonzero value of 
('0a(x)f/'a(x)). The BCS phase has therefore the same 
symmetry as a molecular Bose-Einstein condensate and 
only a crossover occurs. 

The existence of the Ising transition for the case of the 
Bose gas also becomes clear from a consideration of the 
allowed vortices in the two limiting cases. The vortices 
in the state with only a MC are quantized in integer 
multiples of an elementary circulation which is exactly 
one half the quantum of circulation in the state with 
an AC-fMC. Consequently there is a deconfinement-to- 
confinement transition for a pair of elementary vortices 
in the MC with increasing detuning . Experiments in 
the atomic Bose gases thus have the prospect of observ- 
ing this phase transition with an essentially topological 
character, something which has not so far been possible 
in condensed- matter systems. 

The phase diagram that we expect on the basis of the 
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above arguments is shown in Fig. ^ For large posi- 
tive S the molecular energy lies far above the thresh- 
old of the two-atom continuum and the gas consists es- 
sentially completely of atoms, which Bose-Einstein con- 
dense approximately at the ideal gas critical temperature 
T = To = (27rfi2/77i/cB)[n/C(3/2)]2/3, where n is the total 
atomic density of the gas that is assumed to be a con- 
stant throughout the phase diagram. If we lower 6 to- 
wards zero, the number of molecules in the gas increases 
and consequently the critical T decreases monotonously. 
For large negative 5, the molecular energy lies far below 
the threshold of the two-atom continuum, and the gas 
consists solely of molecules that condense at the critical 
temperature Upon increasing 6 towards zero, 

the number of atoms in the gas increases, and the criti- 
cal T for Bose-Einstcin condensation decreases again. Fi- 
nally, our arguments have shown that by increasing 5 at 
a fixed T below the critical T for Bose-Einstein condensa- 
tion of the molecules, a Ising phase transition to a phase 
with both an atomic and a molecular Bose-Einstein con- 
densate occurs. This occurs always at a negative (5, as we 
will see next, when we discuss in detail the quantitative 
determination of the phase diagram. 
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FIG. 1: Phase diagram of a Bose gas with n = 10^ 
a Feshbach resonance as a function of detuning 5 and temper- 
ature T. The solid lines represents the critical T for the tran- 
sitions to the normal state, and the location of the Ising tran- 
sition between the phase with only a molecular Bose-Einstein 
condensate (MC) and the phase with both an atomic and 
a molecular Bose-Einstein condensate (AC-I-MC). The solid 
Unes do not include the effects of the finite lifetime of the 
molecules. The dashed line shows a calculation of the critical 
T for the transitions to the normal state that includes these 
effects. 

Negative detuning. — Up to now we have not made a 
careful distinction between bare and dressed molecules. 
However, for a realistic calculation of the phase diagram 
this distinction is very important. The bare molecule 
corresponds to the bound state of the Feshbach problem 
in the absence of a coupling with the two-atom contin- 



uum, whereas the dressed molecule corresponds to the 
bound state of the full problem including that coupling. 
As a result the energy of the dressed molecule is not 
equal to the detuning, but turns out to be given by 
em ^ (5-h?72(v/l-4Vf?^-l)/2, where rf = g^m^/An^h^ 
and g — hy^2Trahg^BAfi/m determines the coupling be- 
tween the atoms and the bare molecules [2^. The latter 
can be expressed in terms of the background scattering 
length Obg, the magnetic field width AB, and the differ- 
ence between the magnetic moments of a bare molecule 
and two atoms Afi, which are the experimentally know 
parameters of the Feshbach resonance. Furthermore, the 
wave function of the dressed molecule is not equal to the 
wave function of the bare molecule, but contains also am- 
plitude to be in the two-atom continuum. In detail we 

find m 

I Xm; dressed) = ^/Z\Xm; bare) + ^ Ck|k, -k; open), 

k 

(1) 

where the wave-function renormalization factor Z obeys 
l/Z = 1 + r\l1^\trn\- Far from resonance we have that 
Z ~ I and there is essentially no distinction between the 
bare and dressed molecular wave functions. However, 
close to resonance Z becomes very small and the dressed 
molecule actually has almost no amplitude to be in the 
bare molecular state. 

To take the above physics most easily into account 
we introduce creation and annihilation operators for the 
dressed molecules. The grand-canonical hamiltonian for 
the gas then becomes 



H = 



fin - 2/^ 



V'm(a;) 



4m 

2m 2 



V'a(a:) 



dxy^g [V^j„(x)?/>a(a;)V'a(a;) + ilji{x)ipl{x)ipni{x)] 



dx- 



dxT^^ip}^ {x)ipl{x)ipa{x)i(;^{x), 



(2) 



where Tbg — Airaijgh? /m, and Tam and T^ai are the T 
matrices for the scattering of an atom with a dressed 
molecule and for the scattering of two dressed molecules, 
respectively. Near resonance the latter two T matrices 
can actually be expressed in terms of the full atomic scat- 
tering length a{B) = abg[l - AB/{B - Bq)] of the Fes- 
hbach resonance at magnetic field Bq as we show at the 
end of this Letter. 

To find the phase diagram for negative 5, we con- 
sider the phase with only a Bose-Einstein condensate of 
dressed molecules, and perform a quadratic expansion of 
the hamiltonian around the nonzero expectation value 
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('0m (x)) = y'riinc- For atoms with zero momentum, the 
resulting hamiltonian leads to a fluctuation matrix given 
by 



(^^am ^mm)^mc 

4g\/Z^ninc 




mm ) "-mc 



This fluctuation matrix has two positive eigenvalues only 
if 5 is sufficiently negative or, more precisely, if Cm < 
-4gV^V"mc + (STam " ^limOn-mc- For larger 5 an insta- 
bility arises, which physically leads to the formation of an 
atomic condensate in addition to the molecular conden- 
sate that is already present. This condition therefore de- 
termines the position of the quantum (at T = 0) or clas- 
sical (for T > 0) Ising phase transition. The critical T for 
the Bose-Einstein condensation of the dressed molecules 
is of course obtained from the condition = 0. 

Since both conditions are expressed in terms of the 
molecular condensate density rime, we now need to find 
the equation of state of the gas to obtain this condensate 
density as a function of 5 and T . This is complicated 
by the fact that the total atomic density of the gas is 
given by the sum of the density of atoms and of twice 
the density of bare molecules, i.e., not twice the density 
of the dressed molecules. Using the techniques derived in 
Ref.0the calculation can nevertheless be performed and 
we find that in a good approximation the total density 
of bare molecules is 

1 V ^ / ^ ^X'nijii'^mc 1 



Ck + ZTmm'T'mc — 2ftwk \ 



with Ek = ft^k^/2m and fkok = -I- ekT; 



mm ' ^mc 



(3) 



the 



molecular Bogoliubov dispersion. For the atomic density 
we find in a similar manner that 

1 x ^ / 26k ^m (22^am -^mm)^mc 1 



V ^ 



2ek — Cm + (2T'am — Ti-nm)nmc — fli^k 



(4) 



where the dispersion for the atoms obeys hiok = 
V[ek - em/2 + (Tam - T^m/2)n^cV - ^g^Zn^,. Note 
that this dispersion in general has a gap, but becomes 
gapless exactly at the critical condition for the Ising tran- 
sition, as expected. 

Positive detuning. — At positive S no truly stable 
molecular state exists and the previous approach in prin- 
ciple does not apply. That approach was based on the 
fact that for negative 6, the density of states for the bare 
molecules contained a delta function at the energy 
and with strength Z, that corresponded to the dressed 
molecules. For positive 6 the density of states consists. 



however, only of a single broad peak [2l|. Approximating 
that broad peak by a delta function at the energy where 
the density of states has a maximum, we can again use 
the hamiltonian in Eq. but now with 



(5) 



and Z = 1. In this approximation we are neglecting the 
finite lifetime of the bare molecules, that is due to the 
fact that these molecules can now decay into the two- 
atom continuum. At the end of this Letter we, however, 
also discuss a calculation that includes these finite life- 
time effects exactly 19]. The result of that calculation 
is indicated by the dashed line in Fig. 1, which shows 
that the above approximation is indeed rather accurate. 
Since for positive S we only need to determine the crit- 
ical T for Bose-Einstein condensation of the atoms and 
molecules, we here consider only the normal phase of the 
gas. At the critical temperature Tc the atomic density is 
therefore just — 7i(Tc/To)'^/^, whereas the density of 
molecules equals 



= 2V2n f| 



3/2 



53/2 (e 



C(3/2) 



(6) 



where (73/2(2) is the usual Bose function. The desired 
critical T now follows from n = 2n^ + n^. 

Discussion and conclusions. — Up to now we have left 
the molecule-molecule interaction and the atom-molecule 
interaction unspecified. Near resonance they are how- 
ever in first approximation determined by the Feynman 
diagrams given in Fig. |2| Evaluating these diagrams 
at zero external momenta and frequencies, we find that 
7am — 8(7^Z/|em|- Closc to resonance we have that Z = 
AttIi^ /g'^m'^a and |e,„| = h? /ma?. We thus conclude that 
the scattering length for this process is proportional to 
the atom-atom scattering length as Oam = 32a/3. Using a 
similar procedure, the scattering length for the molecule- 
molecule interaction can be shown to be ani,„ = 4a. In 
the case of an atomic Fermi gas we find in the same man- 
ner that Oam = 8a/3 and Omm = ol- This may be com- 
pared to similar results obtained recently by Petrov et al. 
for this case l23l. 



FIG. 2: Feynman diagrams that determine a) the interaction 
between an atom and a dressed molecule and b) the interac- 
tion between two dressed molecules. 
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In the calculation of the phase diagram presented in 
Fig. n we have actually not included the mean-field ef- 
fects due to the atom-molecule and molecule-molecule 
interactions. The reason is that these mean-field effects, 
if we also take into account the unitarity limit of the 
T matrix, are estimated to lead only to relatively small 
shifts. Qualitatively, the most important effect is that 
the Ising transition is, near T = 0, shifted to somewhat 
higher values of 6. In particular, the transition tempera- 
tures for Bose-Einstein condensation are hardly affected. 
As mentioned previously, we have performed a calcula- 
tion of the critical T for Bose-Einstein condensation as 
a function of S, that includes both the effects of the fi- 
nite lifetime of the molecule at positive 6, as well as the 
rogue-dissociation process jl^l for negative S. The results 
are presented by the dashed line in Fig. Note that for 
negative S these effects are nonnegligible only in a small 
region close to resonance, due to the fact that Z becomes 
very small, which implies that there is a rather large spec- 
tral weight in the density of states at positive energies. 
In addition, becomes very small, which implies that 
the positive energy states also get thermally populated. 
For positive S the results are affected more, because of 
the fact that the finite lifetime of the molecules leads to 
a substantial broadening of the density of states, which 
considerably affects the total number of bare molecules 
in the gas. In future work we intend to determine also 
the effects of rogue dissociation on the critical line of the 
Ising transition. We also want to explore the possibility 
that sufficiently close to resonance, the continuous Ising 
transition is actually preempted by a first-order phase 
transition by such effects. 

Field theories for critical fluctuations near the Ising 
transitions at T = and T > can be developed by stan- 
dard methods, using that in the states with (ipni) 7^ the 
Ising order parameter (j> cx i(V'a— V'!)- These field theories 
have the familiar 0^ form, except that at T = there is 
also a marginal coupling to dynamic density fluctuations 
of the superfluid which apparently drives the transition 
weakly first-order [i^ (the T > transition can remain 
second order). The critical precursors of the transition 
lead to a large density of states for collective low energy 
excitations, which will strongly damp single-particle exci- 
tations. We speculate that observation of such enhanced 
sources of dissipation may serve, in addition to direct 
observation of the MC and AC-I-MC, as an experimental 
signature of the Ising transition. 

While the writing of this paper was being completed, 
we became aware of the work of Radzihovsky et al. |26l| , 
who also point out the existence of the Ising phase tran- 
sition, but use the Timmermans model [l^ to determine 
the phase diagram. This work is supported by the Sticht- 
ing voor Fundamenteel Onderzoek der Materie (FOM), 
the Nederlandse Organisatie voor Wetenschaplijk Onder- 



zoek (NWO), and the US National Science Foundation 
grant DMR-0098226. 
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